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ABSTRACT 


Effects  of  horizontal  stochastic  bottom  structure  on  acoustic  intensity  in  an 
isospeed  ocean  have  been  studied  previously  by  the  authors  [J .  Acoust.  Soc. 

Am.  76. »  1445-1455  (1984)].  The  bottom  density  and  sound  speed  were  taken  to 
be  random  and  ray  theory  was  used.  In  this  study  we  add  a  rough  water-bottom 
interface,  consisting  of  large-scale,  two  dimensional  random  facets  upon  which 
small-scale  roughness  may  be  superimposed.  Each  facet  is  assumed  to  possess 
small  random  depth  deviation,  slope,  and  curvature.  Initially,  we  take  acous¬ 
tic  rays  to  be  specularly  reflected  from  the  facet  bottom,  and  derive  expres¬ 
sions  for  the  mean  and  variance  of  incoherent  intensity  at  a  point  receiver 
for  a  transmitted  cw  signal.  The  results  are  sufficiently  general  to  permit 
their  use  with  different  bottom-acoustic  models,  but  we  employ  Mackenzie 
theory  here.  Relative  effects  of  structure  and  topography  are  compared.  Sub¬ 
sequently,  small-scale  roughness  is  superimposed  on  the  facets,  and  the  added 
consequences  of  scattering  are  considered. 


INTRODUCTION 


Acoustical  consequences  of  bottom  characteristics  are  well  known  to  be 
particularly  important  in  shallow  water,  because  of  the  extensive  interaction 
between  sound  waves  and  the  bottom.  This  complication  is  compounded  by  geo¬ 
logical  processes  causing  sedimentation  in  shallow  water,  so  that  there  is 
usually  a  larger  degree  of  lateral  variability  in  bottom  composition  than  in 
the  deep  ocean.1 

Effects  of  bottom  topography  on  acoustic  signals  have  been  studied  ex¬ 
tensively  (e.g.,  Ref.  2).  The  influence  on  sound  propagation  of  horizontal 
ocean  bottom  structure,  which  is  regarded  as  being  characterized  by  its 
density  and  sound  speed,  has  also  received  attention  (e.g..  Ref.  3).  A  few 
investigations  have  discussed  the  situation  when  both  types  of  bottom  varia¬ 
tions  are  assumed  present.4*5  However,  none  of  which  we  are  aware  has  studied 
their  relative  effects,  or  has  used  the  statistical  methods  presented  here. 
Since  structural  and  topographical  properties  are  typically  not  well  known,  it 
is  appropriate  to  model  them  in  a  statistically  random  fashion. 

In  this  paper,  we  use  the  ray  theory  of  propagation,  and  the  vertical 
plane  containing  the  point  source  and  receiver  is  taken  normal  to  a  two-dimen¬ 
sional  rough  bottom.  The  acoustic  field  at  a  point  receiver  is  composed  of 
individual  ray  arrivals  which  are  transmitted  from  an  omnidirectional  cw 
source.  In  much  of  this  study,  ray  paths  between  the  source  and  receiver  are 
determined  by  specular  reflection  from  a  bottom  with  large-scale  roughness  and 
from  a  horizontal  ocean  surface.  Hence,  all  rays  are  coplanar.  Because  we 
are  assuming  an  ocean  of  shallow  depth,  the  water  is  taken  to  have  constant 
sound  speed  and  density.  Our  primary  objective  is  to  determine  statistics  of 
incoherent  intensity  at  a  receiver  in  terms  of  statistics  of  the  bottom  struc¬ 
ture  and  topography. 
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In  an  earlier  paper  by  us,3  the  influence  of  horizontal  structural 
variations  alone  on  acoustic  propagation  was  investigated.  We  use  similar 
methods  here,  and  assume  the  bottom  has  small  random  variations  in  sound  speed 
and  density  in  the  horizontal  direction.  We  generalize  our  previous  investi¬ 
gation  by  permitting  a  rough  water-bottom  interface.  At  first,  we  consider  a 
bottom  composed  of  a  collection  of  large-scale  facets.  Other  studies  have 
used  such  a  facet-bottom  model  to  include  roughness  effects.6*7  Our  facets 
are  characterized  in  part  by  the  locations  of  bottom- reflection  points  of  an 
"ideal"  ray  from  a  horizontal  bottom.  Each  facet  is  assumed  to  have  a  small 
random  inclination  angle  and  depth  deviation.  In  addition,  each  facet  has 
curvature  and  the  radius  of  curvature  is  taken  to  be  large  with  respect  to 
acoustic  wavelength.  Consequently,  the  curvature  of  each  is  modeled  as  a 
random  variable  with  small  magnitude.  Subsequently,  small-scale  roughness  is 
added  to  the  facets. 

In  Sec.  I,  we  present  expressions  for  bottom  loss  and  spreading  loss, 
approximating  both  by  truncated  Maclaurin  series.  In  Sec.  II,  expressions  for 
mean  incoherent  intensity  and  its  variance  are  developed.  Mackenzie  bottom 
loss  and  reasonable  correlation  functions  for  our  bottom  random  variables  are 
incorporated  into  our  equations.  Section  III  presents,  and  provides  discus¬ 
sion  of,  numerical  results  obtained  from  the  model  in  Sec.  II.  In  Sec.  IV,  we 
superimpose  microroughness  onto  our  random  facets,  and  examine  their  conse¬ 
quences  on  intensity.  A  summary  of  principal  results  appears  in  Sec.  Y. 

I.  BASIC  FORMULATION 

Using  the  notation  of  Ref.  3,  we  consider  an  isospeed  sound  channel  of 
mean  depth  H.  The  xy-plane  contains  a  point  source  S  at  (0,hs)  and  a  receiv¬ 
ing  point  R  at  (R,h^).  In  Fig.  1  an  "ideal"  ray  rnj  appears  as  the  dashed 
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curve.  This  ray  has  initial  angle  at  the  source,  the  same  angle  magni- 
tude  at  the  receiver  and  at  each  bottom  reflection.  As  indicated  in  Fig.l, 
all  angles  are  measured  positive  clockwise  from  the  positive  x-direction.  The 
index  j  can  take  on  integer  values  from  1  to  4,  depending  on  the  orientation 
of  the  ray  at  the  source  and  receiver.  In  Fig.  1,  j  has  the  value  one,  which 
corresponds  to  a  ray  with  downward  direction  at  the  source  and  upward 
direction  at  the  receiver.  The  ray  r  .  intercepts  the  mean  bottom  n  times  at 

*  *  J 

the  points  (x:^  ,  H),  i  =  1,  ...,  n. 

Also  shown  in  Fig.  1  is  the  corresponding  "actual"  ray  r  .,  which 

n  J 

appears  as  the  solid  curve.  The  initial  angle  of  this  ray  is  denoted  by  . 
It  intercepts  the  two-dimensional  water-bottom  interface  n  times  at  the  points 
(xj^ ,  i  =  1,  ...,  n.  It  can  be  shown  that,  for  a  bottom  that  deviates 

only  slightly  from  y  =  H,  there  is  an  actual  ray  rnj-  corresponding  to  each 
ideal  ray  r  .  for  each  ordered  pair  (n,j). 

* 1 J 

The  tangent  line  to  the  bottom  at  the  point  (x|j^,  yjj^ )  intercepts  the 
line  x  =  x|^  at  a  depth  that  differs  from  H  by  the  signed  quantity 
Thus,  F.j  is  one  measure  of  the  depth  deviation  of  a  bottom  facet.  It  is 
positive  (or  negative)  when  the  facet  depth  is  greater  (or  less)  than  H  at 

ij 

We  denote  the  small  inclination  angle  of  the  tangent  line  to  the  bottom  at 

(x|j\  y-J^)  by  a!^  and  the  incident  angle  r.  (i.e.,  the  angle  the  incident 

ray  makes  with  the  x-direction  at  the  bottom  reflection)  by  When  the 

^  J 

bottom  facet  in  the  neighborhood  of  the  point  of  reflection  is  concave  down, 
as  is  the  case  shown  in  Fig.  1,  we  take  the  "signed  curvature"  to  be 
positive,  and  when  it  is  concave  up,  is  negative.  Thus,  at  (x(.n^»  yln) ) 
the  bottom  has  a  radius  of  curvature  j =  1/|k^].  Hereafter,  a  primed 
length  quantity  denotes  its  normalization  with  respect  to  the  mean  depth  H 
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(e.g.,  h^  =  h^/H).  We  consider  aj"\  ^  ,  and  =  Hk!^  to  be  random 

variables  of  the  bottom  topography  with  small  magnitudes. 

The  grazing  angle  of  ray  r  .  is  the  angle  between  the  incident  ray 

I  J  II J 

and  the  tangent  line  to  the  bottom  at  (xj^»  y-^),  assuming  specular  reflec- 
tion.  This  angle  is  measured  positive  clockwise  from  a  tangent  line  to  the 
bottom  in  the  ray  direction.  The  grazing  angle  can  be  written  in  terms  of  a 

This  series,  truncated  after 


Maclaurin  series  in  powers  of  a^.n)  and 

TO  1 J 

linear  terms  in  these  variables,  is  given  by 

Jn)  _  ~(n)  (n) 

TO  OJ  M1J  ’ 


(1) 


where 


u(n)  =  2  l  {  — — [2 ( n- £)  +  1  -  Ah'  ]  a^J  + 
ij  4=1  IP  •  R  £J 

•  IJ 

-  2  1i1  a(n)  a(n) 

2  ^  a£0  "  aij  * 


,(n) 


(R 


+  U  .21 

nj  J 


p(n) 1  , 

ij  J 


(2) 


In  Eq.  (2),  u' .  is  the  normalized  version  of  the  total  vertical  distance  u  . 

nJ  nj 

traveled  by  the  mean  ray  r  .,  where 

nj 


“nj  =  2n  H  +  v  hs  -  x  hr 


(3) 


The  symbols  v  and  A  are  parameters  defined  by 

v  =  (-l)j 
and 


(4) 


A  =  (2.5  -  j)/|2.5-j|  . 


(5) 


We  turn  next  to  our  model  for  the  structure  of  the  bottom.  We  consider 

variations  in  the  horizontal  xy-plane,  both  in  the  bottom  density  and 
(n) 

sound  speed  c.  '  at  the  ith  bounce  of  ray  rni .  We  write 


ij 


nj* 

(n).  /  .  'i/i.  (n). 

pi j  /pi  =  (p2/pi J(1  +  eij  )  » 


(6) 
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where  p2  is  the  constant  horizontally  -  averaged  density  of  the  bottom,  pj 
is  the  constant  water  density,  and  is  a  small  random  quantity.  Similarly, 

*  J 


C‘/Cij'  =  tCl/C2)(1  +  ’ 


(7) 


where  Cj  and  c2  are  the  constant  water  sound  speed  and  the  mean  bottom  sound 
speed.  The  quantity  is  another  small  random  variable.  As  discussed  in 


ij 


;(n) 


(n) 


Ref.  3,  it  is  possible  to  approximately  relate  c9  and  <5  .  to  p9  and  e_  , 
respectively,  by  the  equations 


2  ij  ”  2  iJ 


c2(p2)  =  2270.9  -  1194. 4p2  +  474. 6pJ 


and 


(n)  _  _  An) 

eij  =  -  2  6ij 


(8a) 


(8b) 


For  convenience  in  the  expansions  to  follow,  we  assume  that  the  three  topo¬ 
graphical  random  variables,  ,  a!-^  ,  and  ,  and  the  two  structural 

random  variables,  and  6^,  are  of  the  same  small  order  of  magnitude. 

This  enables  us  to  keep  terms  of  like  degree  in  all  variables.  Different 
orders  can,  of  course,  be  assumed,  so  long  as  all  expansions  are  performed  with 
consi stency . 

The  sound  source  S  emits  a  unity  -  amplitude  cw  signal,  which  arrives  at 

the  receiver  R  along  ray  r  .  with  amplitude  a  ..  For  the  large-scale  facet 

n  3  no 

bottom,  the  received  amplitude  is  taken  to  differ  from  unity  because  of  geome¬ 
trical  spreading  loss  and  bottom  loss.  Other  losses  could  be  included,  but 
will  not  be  considered  here.  We  may  write 

a  .  =  A"!/2  B  .  ,  (9) 

nj  no  no 


where  A  .  is  the  spreading  loss  and  B  .  is  the  total  bottom  loss.  In  Ref. 
no  no 

8,  a  spreading  loss  equation  is  derived  for  rays  specularly  reflected  from 
a  bottom  composed  of  large-scale  curved  facets.  Truncated  to  linear  terms. 
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it  is 


Ani  -  H2i(R'2  +7.J)  ♦  2  1  [2^  o‘"»'  ♦  W  ]} 


nj 


nj 


i=l 


ij  10 


(10) 


nj 


,(n) 


where  nlj''  is  given  by 
,(n) 


5J  =  (R*  /“nj  +  0  ["4i  +  2(2n  +  2  “  xhR  "  vhs^'  "  (1  '  vhg) 


(ID 


x  (2n  +  1  -  Ah' )]  . 

R 

This  formula  is  valid  if  n  is  not  too  large.  Note  that,  to  first  order,  the 
slopes  of  the  facets  do  not  contribute  to  A n j .  For  any  particular  ray 


,(n) ' 


are  the  same  for  each  bottom 


r  the  effects  of  the  depth  deviations  V- 
nj  r  ij 

( n ) 1 

reflection  because  the  term  multiplying  IX.  is  independent  of  i.  However, 
the  random  variables  k|j  are  weighted  according  to  which  bottom  is  consid¬ 
ered;  the  curvature  of  bottom  facets  near  midrange  has  the  greatest  effect, 
and  that  near  the  source  and  receiver  has  the  least  influence. 

In  Eq.  (9),  the  decrease  B  .  in  amplitude  due  to  all  the  bottom  losses 

■  *  J 


of  ray  rnj  is 


where  is  the  reflection  coefficient  of  r  .  at  its  i^1  bottom  bounce.  We 

J  (n)  ° 

assume  that  the  coefficient  B[..  depends  on  both  grazing  angle  and  bottom 

*  si 

( n ) 

structural  properties.  It  is  then  possible  to  expand  B. .  in  a  Maclaurin 

'  J 

series  with  respect  to  the  random  variations  of  these  quantities.  To  first- 
order  terms,  there  is  then  an  expression  with  the  form 

(13) 


B  . 
nj 


",  e<"> 

nj 


(12) 


i=l 


B<n)  =  A  •  +  B  ,t!n)  ♦  C  -6ln)  ♦  D 
ij  nj  nj  ij  nj  nj  njnj 


(n)  . 


where  u. .  is  given  by  Eq.  (2).  The  coefficients  of  the  random  variables  in 

*  J 

Eq.  (13)  depend  on  the  reflection  theory  used  and  the  mean  bottom  values 
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selected. 


II.  INTENSITY  MOMENTS 


We  consider  two  measures  of  the  sound  field  at  the  receiver.  These  are 

the  mean  and  variance  of  the  incoherent  intensity  I,  given  by 

N  4  o 

I  =  l  la.  (14) 

n=0  j=l 

Because  the  losses  experienced  by  successive  ray  arrivals  grow  rapidly,  we 
need  consider  only  those  rays  experiencing  some  number  N  or  fewer  bottom  re¬ 
flections.  Specifically,  we  elect  to  neglect  ray  arrivals  which  have  ampli¬ 
tudes  less  than  1%  of  the  arrival  with  largest  amplitude. 

We  wish  to  determine  expressions  for  the  mean  and  variance  of  I, 

u(I)  =  E(I)  (15) 

and 


a2(I)  =  E(I2)  -  E2(I)  , 


(16) 


where  E  denotes  expectation  and  I  is  a  function  of  the  five  types  of  random 
variables  introduced  in  Sec.  I.  We  assume  that  the  mean  of  each  type  is  zero, 

.  (17) 


and  that  each  type  has  a  constant  standard  deviation, 


°hj  )  =  °e  » 

(18a) 

)  =  o. , 
v  1J  ;  6  * 

(18b) 

rdn) 1  ^ 

) =  V 

(18c) 

(An), 

o[a..  )  =  oa  , 

( 18d) 

and 

r  (n)\ 

°hj  3  -  v  • 

(18e) 

-  8  - 


We  also  assume  that  any  two  different  topographical  variables  are  uncorrelated, 
and  that  topographical  and  structural  random  variables  are  not  correlated.  For 
example,  the  concavity  at  a  bottom  reflection  is  taken  to  be  uncorrelated  with 
both  the  depth  deviation  and  the  density  there.  However,  we  do  allow  each  ran¬ 
dom  variable  to  be  correlated  with  itself  at  two  separate  bottom  reflections. 
Thus, 

Efe1 

10 


FrE(n)  (m)  'i  _  2  ,  (n)  (mK 

Eleij  e*k  J  -  °e  Celeij  •  e£k  J  * 


r  c  x ( n )  Am)  >  2  ( Jn)  Am)  > 
El6ij  Sik  >  -  °&  C6l*ij  '  S*k  > 


E(P<">'  p<™>' 
'■10 


'i  _  2  o  rf)(n)'  fl(m) 

j  '  Ak  J  "  °V'  V'  ^  ij  *  P£k 


)  > 


rf.tn)  (m) x  _  2  r  (  (n)  (m) , 

Elaij  a*k  J  =  °a  Va1j  •  akk  )  ■ 


and 


r-  f  \  >' 

Ehj 


(n) *  (m) ’ 


'£k 


)-  -c  , 

K  K  1 J 


(n) ’  (m) ’ 

c£k 


)  • 


(19a) 

(19b) 

(19c) 
( 19d) 

(19e) 


In  Eqs.  (19),  any  function  C(a,b)  represents  the  correlation  coefficient  of  the 

random  quantities  a  and  b.  Also,  the  structural  random  variables  and  6^ 

10  io 

are  correlated  as  in  Ref.  3,  so  that 


cr  'n 

Eh.i 


(n)  Am) 


0  £k 


1  _  _  _  r  <•  (n)  (m)  \ 

1  Cc6^ei.i  »  6£k  j  • 


->  o  .  I  e .  • 

e  6  1J 

Using  the  particular  relationship  given  by  Eq.  (8b),  we  may  write 

2  _  f  An)  _(m)  >  _  .  2 


and 


2  (  An)  (m)  x  2  ,  (n)  (m)  , 

°6  C6(6ij  >  6tk  )  =  1/4  °e  CeCeij  •  £*k  )  * 


rr  a  C  f 1  -  1/9  ^  C  f 1 

°e°6  e6^ei j  ’  6£k  J  “  -1/2  °£  M£ij  ’  e£k  )  ’ 


(n)  (m) 


( 1 9  f ) 

(20a) 

(20b) 


so  that  we  have  only  one  correlation  function  C  for  the  structure. 

e 

It  is  now  possible  to  determine  w( I )  and  o2(I).  Keeping  through  first- 
order  terms  in  the  random  variables  for  the  mean,  we  have 
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y( I )  =  E  E  A^n/[H2(R'2  +  u’  2)]  .  (21) 

n=0  j=l  J  J 

In  Eq.  (21),  A  .  and  H2 (R ' 2  +  u'  z)  represent  the  leading-order  terms  of  the 


r.2 

nj  v  nj 

bottom  loss  and  spreading  loss,  respectively  (see  Eqs.  (13)  and  (10)).  No 


first-order  terms  appear  in  Eq.  (21)  due  to  the  assumptions  of  Eq.  (17).  We 
keep  through  second-order  terms  in  the  variance  so  that  the  standard  deviation 
is  correct  to  first  order.  Therefore,  for  the  variance,  we  obtain  an  expres¬ 
sion  of  the  form 


cT(I)  =  V  o2  +  ,  a2,  +  V  a2  +  V 

e  e  V  v  a  a 


2 

■  o  , 
<  < 


(22) 


The  quantities  V£o^  ,  Vp,c^,  ,  ,  and  V^,a2,  in  Eq.  (22)  represent  symboli¬ 

cally  the  contributions  to  the  variance  from  the  bottom  structure,  depth  devia¬ 
tion,  slope,  and  curvature.  After  much  manipulation,  the  coefficients  can  be 
found  as  the  following  formulas: 


N  N  4  4 

V  =  !_  E  l  E  E 


.2n-l  .2m-l 
Anj  Amk 


H 


*+  n= 


_  fB  -B  ,  -  B  -C 

I-  m  ml-  m  r 


n=l  m=l  j=l  k=l  (R'2  +  u;?)(R’2  +  ^2)  nj  mk  nj  mk 


n  m 


+  I  C  .C  .  ]  E  EC  (e(n) ,  e(m)  ) 
4  nj  mk  ^  ^  G  lj  tk 


(23a) 


N  N  4  4 
V  .=  1°  E  E  E  E 
V 


»2n  .2 
'nj 


*2 


H4  n=l  m=l  j=l  k=l  (R'2  +  u;?)2(R’2  +  u’?)2 

‘  *  J 


rR'Dnj  n  _  1 

1  A^r 


nj 


'nj 


x  fR  Dmk  m  _  1  y  ™  r  ro(n)',o(m)\ 

Her  v)  in  v“  1  ’ 


(23b) 


N  N  4  4 

V  =  £_  E  E  E  E 


•  2n-l  .2m-l  n 
Anj  ^mk  Dnj 


H*  "=1  -1  3=1  k=1  “nj  “;k(R'2  +  “nj)(R'Z  +  i=1  1=1 


n  m 

E  E  { [2( n-i ) 


+  l]vh;  +  (21-l)xh;H[2(»-t)  +  l]vn;  +  (2k-l)xh;l  a4!1)  , 


(23c) 
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and 


N  N  4  4 

V  ,=  V  z  z  z  z 


.2n  .2m 
^nj  ^mk 


h*  n=1  npl  j=1  k=1  Kj  uri JR'2  +  u^kr'2  +  ^k2) 


(23d) 


n  m 
z  z 
i=l  1=1 


,(n)  (m) 

ij  n£k 


r  r  (n) '  (m) '  ^ 

x  £  £  a--  0..  CK.  [Klj  Kjk  )  . 


Note  that  the  variance  depends  only  on  the  first-order  coefficients  B  ..  C  .. 

nj  nj 

and  D  .  of  the  random  variables  in  Eq.  (13),  and  not  on  any  second-order 

1 1 J 

coefficients.  That  is,  even  if  the  coefficients  of  second-order  terms  were 
somehow  known  in  Eq.  (13),  these  would  not  appear  in  Eqs.  (22)  and  (23),  due 
to  the  averaging  process. 

Investigation  of  Eqs.  (21)  and  (22)  for  the  mean  and  variance  of  inten¬ 
sity  requires  specification  of  a  bottom  reflection  model.  We  choose  here  to 
use  Mackenzie's  model,9  in  which  the  reflection  coefficient  is  given  by 


( n)  r ,,  .  (n) ,,  l,,i/2 

B-j j  =  [(h+  -  a  sin  )2  +  h2] 


(24a) 


(n) 


*  [ ( h+  +  o  sin  )2  +  hf  ] 


1-1/2 


where 


h+  =  {±  g  +  [(a/3)2  +  g2]1/2}172  , 


9  =  i  [1  -  (Vs)2  cos2  '('ij'  -  (a/e)2]  , 


(n) 


(24b) 

(24c) 


and 


(n)  . 

0=p(j  Vj/PjC,  . 


(24d) 


The  quantity  a  is  the  impedance  and  v,  is  the  phase  speed  of  acoustic  waves  in 
the  bottom  and,  to  an  excellent  approximation,  is  given  by  v2  =  •  The 

quantity  a/ 3  is  the  ratio  of  the  imaginary  to  real  parts  of  the  wave  number  of 
acoustic  waves  in  the  bottom.  This  term  provides  a  dissipative  mechanism  for 
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sound  which  enters  the  bottom.  In  Ref.  3  it  is  explained  how  values  for  a/e 
are  determined  from  formulas  derived  by  Hamilton.10 

Using  Mackenzie  theory,  we  obtain  expressions  for  the  coefficients  Anj , 

B  C  and  D  •  which  appear  in  Eq.  (13).  They  are  calculated  to  be 
nj *  nj ’  ni  rr  ^ 


and 


where 


nj 


A  •  =  M./M  •  , 
nj  nj  nj 


B  .  *  Q  .  (S  .  -  P  .  )/[(M+.  )3M"  ]  , 
m  xm  ^  m  nj  }  L  v  nj  '  n J  9 


'nj  nj  v  nj 
“2 


(25a) 

(25b) 


=  Q  .  c2  (P  .  +  S  .  -  4T  . )  +  Q  .  (S  .  -  P  .) 

yni  n.i  n.i  nj 1  ynj  *■  no  nj ; 


nj  nj 

x  [p  VoVor1  cos2 "e^  , 
L  m  m  >  v  nj  '  J  oj  * 


nj  ^nj 

^  .  (K  . 

nj  nj 


(25c) 


D  .  =  R'Q  .fS  .(1  -  1/p2)  -  P  .  +  (S2./P2  +  4T  .c  sin|e^l)/P  .1 
nj  ynj L  nj k  J  nj  k  nj  nj  nj  >'  nj J 


x  [u'.(M  .)3M'.]_1  , 
L  nj v  nj  ;  nj  J 


nj 
(25d) 


M".  *  [P  .  ±  Q  .  +  S  .1 
nj  L  nj  ynj  nj J 

P  .  =  2  fa2/ e2  +  T  .2) 
nj  v  m  J 


1/2 

1/2 


nj 


ITnj  +  («2^2  +  T„/)*/2]l/2  *fn|e*3»| 


S  .  =  7  sin2  el">  , 
nj  oj  * 

Tnj  =  i  t1  -  c2  cosl^  -  aW]  , 

o  =  c  p  , 

C  =  c2/c1  , 

■p  =  p2/px  • 


(25e) 
( 25  f ) 

(25g) 

{ 2  5  h ) 

(251 ) 

(25j ) 
(25k) 


and 


(25£) 
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The  above  formulas  reduce  to  those  in  Ref.  3  where  the  bottom  was  horizontal. 
The  notation  has  been  changed  slightly,  however. 

To  evaluate  the  intensity  variance  in  Eq.  (22),  we  also  need  to  specify 
correlation  coefficients  for  the  random  variables.  We  choose  the  Gaussian 
forms 


(n) 

ij  ’ 

( m )  -v  r  r  ( n)  (m)  \2  „  2 1 

e£k  )  exp  ^'(xij  "  x£k  )  /LJ  » 

(26a) 

(n)  > 

ij 

=  pvn  r_ry(n)  ( m )  >. 2  ^  •, 

*  £k  '  p  *■  ^xi  j  ”  X£k  ^  i  * 

(26b) 

(n) 
ij  * 

)  *  exP  I-(Nj  -  \k  )  /Lal  ' 

(26c) 

(n) 

K.  . 

1  J 

'  K(m)'  1  =  exp  r-fx(n)  x(m))2/i  2  1 

’  K£k  }  6XP  l  ^ij  X£k  )  /l<'  J  » 

(26d) 

and 


Here,  |x|n^  -  |  is  the  horizontal  distance  between  the  ith  bottom  reflection 

l  J 

of  ray  r  .  and  the  i  th  bottom  bounce  of  ray  r 
nj  J 


k.  The  parameters  L  ,  L  L  , 

t  V  a 


L^,  are  the  correlation  lengths  associated  with  each  random  variable,  and  each 
represents  the  value  of  |x*"}  -  x^}|  where  c  takes  the  value  1/e.  The  correl¬ 
ation  function  for  the  bottom  sound  speed  variation  6^  is  related  to  C  bv 
Eqs.  (20).  Of  course,  other  forms  for  the  correlation  functions  could  be  used, 
if  desired  or  if  suggested  by  bottom  observations. 


III.  NUMERICAL  RESULTS 


We  now  investigate  the  behavior  of  the  mean  and  variance  of  the  received 
incoherent  intensity,  Eqs.  (21)  and  (22),  respectively.  We  incorporate  expres¬ 
sions  from  Eqs.  (23)  into  the  variance  and  use  the  Mackenzie  bottom  model 
formulas,  Eqs.  (24)  and  (25),  for  the  coefficients  of  bottom  loss  and  spreading 
loss.  For  the  correlation  coefficients  in  the  variance,  we  use  the  Gaussian 
forms  shown  in  Eqs.  (26).  As  indicated  in  the  theoretical  development,  our 
expressions  apply  to  relatively  large-scale  topographical  variations. 
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Graphs  of  the  moment  ratio 

Mr  =  10  log1Q  p(I)/o(l)  (27) 

will  appear  in  our  figures.  This  is,  in  effect,  a  signal-to-noise  ratio  in  dB. 
Small  (i.e.,  large  negative)  values  for  MR  mean  that  the  fluctuation  in  inten¬ 
sity  is  large  relative  to  its  mean.  Conversely,  large  (positive)  values 
correspond  to  the  standard  deviation  being  small  compared  to  the  mean,  i.e.  a 
stronger  "signal." 

For  convenience  in  our  numerical  calculations,  we  place  the  source  and 
receiver  on  the  ocean  surface.  Of  course,  other  configurations  could  be  used 
without  difficulty,  if  desired.  The  density  and  sound- speed  values  in  the 
water  are  taken  to  be  =  1.025  g/cm3  and  ^  =  1523  m/s.  We  found,  for  the 

parameter  values  chosen,  that  a  maximum  number  N  of  at  most  7  bottom  reflec¬ 
tions  needed  to  be  included  for  any  aspect  ratio  R' ,  according  to  the  criterion 
presented  in  Sec.  II.  Rays  with  a  greater  number  of  bottom  reflections  were 
shown  to  negligibly  affect  the  results.  From  Eq.  (23c),  it  is  apparent  that 

for  source  and  receiver  on  the  surface  (h^  =  hj.  =  0) ,  the  contribution  from  the 
bottom  slope  vanishes  because  V,  =  0.  From  numerical  calculations  for  other 

a 

source  and  receiver  locations,  we  found  that  the  contributions  from  the  bottom 
slope  terms  is  always  insignificant  compared  to  the  structural  and  curvature 

contributions.  The  quantity  can  also  be  demonstrated  to  be  relatively 
small  compared  to  V£  and  V^,. 

Figure  2  shows  moment  ratio,  Eq.  (27),  versus  mean  bottom  density  p  . 

The  values  used  for  p2  are  typical  of  bottoms  on  the  continental  shelf.  Based 
on  data  from  various  sources,11  we  have  chosen  the  values  of  the  standard  de¬ 
viations  for  our  random  variables  to  be  0.1  (for  0,0^,,  and  o  )  and  0.2 

6  V  1c1 

(for  a  )•  These  are  consistent  with  our  earlier  assumption  that  the  random 
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variables  are  of  the  same  order  of  magnitude  and  with  Eq.  (8b).  We  have  kept 
the  aspect  ratio  constant  (R‘  =20)  in  this  figure,  and  each  curve  corresponds 
to  different  values  of  the  two  normalized  (with  respect  to  H)  correlation 
lengths  and  L^, .  Also,  we  have  taken  the  correlation  lengths  of  the  two 
topographical  random  variables,  L  ,  and  L  , ,  to  be  equal. 

{/  1C 

In  the  figure,  and  L^, ,  take  on  values  of  zero  and  infinity.  When 
L\  is  near  zero,  it  means  there  is  little  correlation  between  bottom  topo¬ 
graphy  at  any  two  given  points.  On  the  other  hand,  when  L\  is  very  large, 
the  topography  is  correlated  even  at  widely  separated  bottom  locations. 

Similar  statements  apply  for  the  normalized  correlation  length  of  bottom 
structure,  L^.  It  is  apparent  from  Fig.  2  that  varying  the  normalized  topo¬ 
graphical  correlation  length  L1 ,  has  a  much  greater  effect  on  MD  than  varying 
L' •  For  example,  at  p2  =  2.1  g/cm3,  when  L\  increases  from  zero  to  infinity, 
there  is  about  a  4  dB  change  in  MD.  In  contrast,  varying  L1  has  very  little 
effect.  The  quantity  L'  has  some  influence  only  when  the  mean  bottom  density 
is  small.  Indeed,  this  is  where  the  Mackenzie  reflection  coefficient  tends  to 
vary  most  rapidly  with  changes  in  bottom  density  or  sound  speed.  However,  at 
no  p0  value  does  varying  L'  cause  more  than  a  1  dB  change  in  MD.  This  figure 
suggests  that  the  moment  ratio  is  insensitive  to  the  precise  value  of  the  hor¬ 
izontal  correlation  length  for  density  or  sound-speed  variations. 

When  the  normalized  correlation  length  L',in  Eq.  (26d)  increases,  the 

K 

correlation  coefficient  C  ,  increases  as  well.  Thus,  for  small  L ' , ,  the 

K  K 

correlation  between  the  bottom  concavity  at  two  points  is  small,  and  the 
fluctuations  in  concavity  at  these  points  tend  to  cancel.  Consequently,  a(I) 
decreases,  causing  the  moment  ratio  to  increase  as  L' ,  decreases.  For  larger 
L^,,  there  is  more  of  a  relationship  between  the  concavities  at  two  bottom 
points.  The  reduction  in  fluctuations  will  not  be  as  great  as  for  small  L'i* 
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Hence,  a(I)  grows  as  L1 ,  increases  and,  therefore,  M  decreases.  The  same 

k  K 

type  of  argument  applies  for  any  of  the  other  random  variables  and  their 
corresponding  correlation  lengths.  However,  as  suggested  by  Fig.  2,  it  is  the 
variations  of  L\  which  have  the  largest  effects  on  our  model.  With  this 
picture  of  the  effects  of  different  values  for  correlation  lengths,  we  proceed 
to  investigate  the  influence  of  other  parameters.  In  all  subsequent  figures 
of  this  paper,  we  will  take  all  correlation  lengths  as  approximately  zero. 

Figure  3  illustrates  the  behavior  of  the  moment  ratio  versus  mean  bottom 
density  for  three  values  of  the  aspect  ratio  R' .  For  a  given  value  of  R' ,  we 
observe  a  general  upward  trend  in  M  as  p  increases.  This  same  trend  can  be 

K  ^ 

seen  in  Fig.  2,  of  course.  For  low-density  bottoms,  the  Mackenzie  reflection 
coefficient  tends  to  fluctuate  more  as  p2  varies  than  for  larger  p^.  This 
causes  a(I)  to  be  large,  so  that  M  in  dB  becomes  more  negative.  For  example, 
when  R'  =5,  there  is  approximately  8  dB  difference  between  the  values  of  MR 
at  p2  =  1.5  g/cm3  and  p2  =  2.1  g/cm3.  It  might  be  anticipated  that  the  var¬ 
iance  of  the  intensity  would  become  smaller  as  R'  decreased  and  thus  that  M 

R 

would  increase.  This  is  true  for  high-density  bottoms,  but  not  for  lower 
values  of  p^.  At  p2  =  1.6  g/cm3,  for  example,  MR  is  larger  for  R1  =  10  than 
for  R1  =  5.  As  was  stated  in  Ref.  3,  low-density  bottoms  have  a  less  well- 
defined  behavior  as  compared  to  higher  density  ones. 

In  Fig.  4,  we  show  the  relative  effects  on  moment  ratio  of  random  bottom 
structure,  random  bottom  topography,  and  a  combination  of  both  features.  The 
quantity  MR  is  plotted  against  the  aspect  ratio,  with  p2  held  constant  at 
1.8  g/cm3.  The  heavy  solid  curve  represents  the  combined  effects  of  topography 
and  structure  on  moment  ratio  ( i . e . ,  all  standard  deviations  in  Eq.  (22)  are 
non-zero).  The  dashed  curve  corresponds  to  a  horizontal  bottom  with  random 
structure  (i.e.,  ap,=  =  o^,  =  0  in  Eq.  (22)).  Similarly,  the  light  solid 
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curve  corresponds  to  a  structurally  homogeneous  bottom  with  random  topography 
(i.e.,  =  0  in  Eq.  (22)).  The  relatively  large  (i.e.,  small  negative)  values 

^  on  the  heavy  solid  curve  for  small  R'  indicate  relatively  small  intensity 
fluctuations.  On  the  other  hand,  the  small  (i.e.,  large  negative)  values  of 
Mr  for  large  R'  correspond  to  large  variations  of  the  intensity  about  its 
mean.  This  same  trend  will  be  seen  in  Figs.  5  and  7. 

It  is  also  apparent  in  Fig.  4  that  for  small  aspect  ratios  (approximately 

R'  <  5) ,  the  random  bottom  structure  makes  a  greater  contribution  to  M  than 

R 

does  the  random  bottom  roughness.  For  larger  ratios,  however,  the  opposite  is 
true.  The  dashed  curve,  corresponding  to  a  horizontal  bottom  with  random 
structure  only,  tends  to  oscillate.  This  is  due  to  the  large  variability  in 
the  bottom  reflection  coefficient  with  grazing  angle.  The  light  solid  curve, 
which  corresponds  to  a  structurally  homogeneous  bottom  with  roughness,  is  dom¬ 
inated  by  curvature  effects.  Because  the  bottom  reflection  coefficient  does 
not  depend  on  the  interface  curvature,  this  curve  monotonically  decreases  with 
increasing  R' .  The  bottom  curvature  tends  to  focus  and  defocus  ray  bundles 
through  the  spreading  loss  contribution.  This  causes  ct(I)  to  be  greater,  and 
the  effect  increases  as  the  grazing  angle  decreases,  i.e.,  as  R'  increases. 

Figure  5  is  a  graph  which  can  be  compared  with  Fig.  4.  The  mean  bottom 
density  is  decreased  from  p2  =  1.8  g/cm3  (e.g.,  fine  sand)  in  Fig.  4  to  p2  = 

1.6  g/cm3  (e.g.,  silty  clay)  in  Fig.  5.  The  general  shapes  of  the  corresponding 
curves  in  each  figure  are  broadly  similar.  However,  in  the  curves  of  Fig.  5, 
the  moment  ratio  becomes  more  negative.  This  is  due  to  the  influence  of  the 
random  structure,  not  to  that  of  the  topography.  As  we  discussed  in  connection 
with  Fig.  3,  a  lower-density  bottom  usually  results  in  a  larger  value  for  o(I). 
The  smaller  mean  density  value  used  for  Fig.  5  also  has  the  effect  of  lengthen¬ 
ing  the  interval  on  the  R1  axis  in  which  structure  contributes  more  to  MD  than 
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topography.  This  interval  extends  to  about  R'  =  9  in  Fig.  5,  compared  with 
about  R '  =5  in  Fig.  4. 

A  three-dimensional  graph  is  shown  in  Fig.  6,  in  which  we  have  varied 

the  standard  deviations  of  the  two  dominant  random  variables,  o  and  a  , . 

e  k 

This  perspective  figure  shows  how  these  quantities  influence  the  moment  ratio. 
The  aspect  ratio  and  mean  bottom  density  are  held  constant  at  10  and  1.8  g/cm3 
Each  of  the  standard  deviations  is  varied  from  0  to  0.2.  When  a  ,  is  near 

K 

zero,  increasing  o  has  the  effect  of  decreasing  MD.  As  o  ,  becomes  closer  to 

0.2  however,  varying  o  has  little  effect  on  Mn.  For  a  constant  value  of  o  , 

increasing  o  ,  always  causes  a  decrease  in  MD.  These  types  of  behavior  occur 
because  at  R1  =  10,  the  topography  of  the  bottom  has  a  larger  influence  on  the 
intensity  variance  than  does  the  structure.  If  R'  were  chosen,  for  example, 
to  be  less  than  five,  then  varying  would  have  a  greater  effect  than 
changing  o  , ,  because  then  structure  dominates  over  topography.  This  was  seen 
previously  in  Fig.  4. 

Figure  7  illustrates  how  the  aspect  ratio  R'  and  mean  bottom  density  p2 
affect  moment  ratio.  In  fact.  Figs.  3,  4,  and  5  are  cross-sections  of  Fig.  7, 
but  were  introduced  first  because  of  their  greater  ease  of  interpretation. 

For  a  given  density,  the  dominant  n  =  1  (one  reflection)  ray  has  a  specific  R' 
value,  or  grazing-angle  value,  at  which  the  Mackenzie  reflection  coefficient 
drops  from  a  value  near  one  to  a  much  lower  value.  This  occurs  near  the 
critical  angle  of  Rayleigh  reflection  theory.12  In  turn,  a  large  a(I)  value  is 
produced,  and  hence,  a  large  negative  value  for  M  .  The  indentation  in  the 

K 

surface  of  Fig.  7  is  the  general  consequence  of  this  feature.  As  p2  increases 
the  value  of  R'  at  which  this  effect  occurs  decreases,  since  the  grazing 
angle  increases.  For  rays  with  n  >  1,  the  same  effect  occurs,  but  these  rays 
are  attenuated  by  each  bottom  reflection  and,  consequently,  do  not  have  so 
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large  an  influence  on  received  intensity.  It  can  also  be  seen  from  Fig.  7 
that,  for  any  p2value,  the  moment  ratio  experiences  an  overall  decrease  as  R' 
increases  form  1  to  20.  This  decrease  is  most  dramatic  for  lower  density 
bottoms. 

IV.  SMALL-SCALE  ROUGHNESS 

In  previous  sections,  we  assumed  specular  reflection  from  a  facet  bottom, 
implying  facet  lengths  which  are  large  in  acoustic  wave-wengths.  In  this  sec¬ 
tion,  we  enhance  our  bottom  model  by  superimposing  small-scale  roughness  on  the 
facets.  It  is  not  uncommon  for  ocean  bottoms  to  be  modeled  using  two  degrees 
of  roughness. 13-16 

Eckart17  gave  a  Kirchhoff  solution  for  the  acoustic  signal  reflected  by 
a  rough  surface  which  was  later  modified  by  Clay.18  We  use  Clay's  method  here 
to  describe  the  sound  scattered  by  the  small-scale  roughness.  A  more  general 
approach,  which  reduces  to  Clay's  results,  is  given  by  Kuperman.19 

We  assume  that  the  dimensions  of  the  scattering  area  are  much  greater 

than  the  acoustic  wavelength  X,  and  these  dimensions  are  in  turn  much  less 

than  the  distance  to  the  source  or  to  the  receiver.  Let  y  =  ci^(x)  be  the 

i  J 

deviation  of  the  irregular  surface  from  the  smooth  facet  which  reflects  ray  r^ 
at  its  i^  bounce.  We  impose  the  usual  conditions  on  this  function: 

U-jVxj  «  1  (28a) 

and 

|<u!^/dx|  «  l  .  (28b) 

Equation  (28a)  states  that  the  magnitude  of  the  small-scale  roughness  is  small 
compared  with  the  acoustic  wavelength.  The  assumption  given  in  Eq.  (28b)  is 
that  the  slope  of  y  =  ?j^(x)  is  never  large.  Because  of  the  restriction  in 
Eq.  (28b),  the  reflection  coefficient  on  any  given  facet  is  nearly  constant, 
with  respect  to  grazing  angle,  for  reflections  in  the  specular  direction. 
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Moreover,  it  has  been  shown  that  curvature  of  the  small-scale  roughness  does 
not  significantly  affect  the  scattering  of  the  acoustic  wave,20  so  we  ignore 
this  feature  here.  If  we  assume  to  be  normally  distributed  with  mean  zero 

over  the  area  insonified  by  a  ray  bundle,  then  the  coherent  reflection  coeffi¬ 
cient,  i.e.,  the  reflection  coefficient  in  the  specular  direction,  is  9iven 
by18 


'(n) 


B(n) 


exp 


[-  8*n2a^i 


S1IV 


♦<?>] 
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(29) 


Here,  b!?  is  the  reflection  coefficient  given  in  Eq.  (13),  and  is  the 

grazing  angle,  defined  prior  to  Eq.  (1),  that  the  incident  ray  r  makes  with 

nj 

the  smooth  facet  at  its  ith  bottom  reflection.  The  quantity  c  is  the  root- 

£ 1 

mean-square  displacement  of  the  surface  from  the  smooth  facet  upon  which 

it  is  superimposed,  divided  by  the  acoustic  wavelength  X.  We  assume  a  (  is 

small  by  using  Eq.  (28a)  and  approximate  8  —  ^  by  ? 

1  J 


o(n)  n(^)  r*  o  ?  ?  *  ?  (n)\ 

B.  .  =  8.  .  (1  -  8.2  Sl»^  )  , 


(30) 


to  first-order  terms  in  o2, . 

c 

In  order  to  incorporate  Eq.  (30)  into  our  previous  development,  we  assume 
that  scattered  waves  at  the  receiver  vanish  because  of  their  random  phases. 

Consequently,  only  the  wave  reflected  in  the  specular  direction  remains.  Prior 
to  inserting  B  into  our  equations  for  the  mean  and  variance  of  received 

intensity,  we  regard  o^,  to  be  the  same  on  each  facet.  Also,  we  assume  they 
are  of  the  same  magnitude  as  a^.  The  mean  intensity  then  becomes 


E(I) 
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c  nj  . 

A2  (R'2  +  u'2) 
nj  v  nj  J  -* 


(31) 


Because  of  our  assumptions  on  o^( ,  the  intensity  variance  remains  unchanged 
from  Eq.  (22),  at  least  to  second-order  terms  in  our  random  variables. 
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In  Fig.  8  we  plot  moment  ratio  Mp  versus  aspect  ratio  R'  for  three  values 
of  ,  0,  0.1,  and  0.2.  The  bottom  has  a  mean  density  of  1.8  g/cm3.  The  a^,=  ( 
curve  is  the  same  as  the  heavy  solid  curve  in  Fig.  4,  since  this  value  of  o  , 


corresponds  to  facets  with  no  superimposed  small-scale  roughness.  For  R'  =  1, 

varying  o  ,  from  zero  to  0.2  produces  about  a  3  dB  decrease  in  M  .  However,  the 
C  R 

small-scale  roughness  in  this  model  is  essentially  inconsequential  for  R'  greater 
than  about  four.  In  other  words,  when  the  grazing  angle  of  the  dominant  rays 
becomes  relatively  small,  the  small  surface  variation  q1?  has  little  effect,  as 
is  well  known.18  An  increase  in  o  ,  causes  a  decrease  in  MD  for  every  value 
of  R'  in  Fig.  8.  This  means  that  there  is  a  tendency  toward  increased  fluctua¬ 
tion  of  intensity  about  its  mean  value  as  a^,  increases. 


V .  SUMMARY 

This  paper  is  concerned  with  the  relative  effects  of  bottom  structure  and 
topography  on  the  intensity  of  a  received  acoustic  signal.  We  use  ray  theory  in 
a  shallow-water  duct  having  constant  sound  speed  and  density.  Both  the  bottom 
sound  speed  and  density  are  taken  to  be  random  in  the  horizontal  direction.  The 
water-bottom  interface  consists  initially  of  an  ensemble  of  large-scale  two 
dimensional  facets.  Each  facet  is  assumed  to  possess  a  small  random  depth  devia¬ 
tion  about  a  mean  horizontal  bottom,  as  well  as  small  random  slope  and  curvature. 
Acoustic  rays  are  specularly  reflected  from  the  facet  bottom. 

Using  perturbation  expansions,  we  develop  expressions  for  bottom  loss  and 
spreading  loss.  These  results  are  used  to  derive  formulas  for  the  mean  and  vari¬ 
ance  of  the  incoherent  intensity  at  a  point  receiver  for  a  transmitted  cw  signal. 
Stochastic  averaging  over  the  bottom  randomness  is  performed  to  obtain  the  mean 
and  variance  of  intensity.  The  results  are  sufficiently  general  to  permit  their 
use  with  different  bottom- acoustic  models.  The  intensity  moments  are  modeled  to 
contain  arbitrary  correlation  coefficients  for  each  of  the  random  variables.  To 


-  21  - 


illustrate  our  results,  we  choose  to  use  a  Mackenzie  bottom  model  and  correl¬ 
ation  coefficients  of  Gaussian  form. 

We  are  able  to  draw  a  number  of  conclusions  from  our  results.  For 
example,  we  show  that  the  standard  deviation  of  intensity  increases  as  the 
correlation  length  of  bottom  structure  or  topography  increases.  This  results 
from  a  cancelling  of  bottom  randomness  effects.  Another  conclusion  is  that 
varying  the  correlation  length  of  the  topography  affects  the  moment  ratio, 
the  dB  ratio  of  mean  to  standard  deviation  of  intensity,  more  than  varying 
that  of  the  structure.  Moreover,  varying  the  correlation  length  of  any  of 
the  random  bottom  quantities  does  not  affect  the  moment  ratio  as  much  as  does 
varying  the  aspect  (range  to  depth)  ratio  R'  does.  Also,  in  general,  for  a 
given  R‘  value,  the  moment  ratio  increases  as  the  mean  bottom  density  increases 
When  considering  with  high  density  bottoms,  the  moment  ratio  tends  to  decrease 
as  the  aspect  ratio  increases.  However,  this  same  conclusion  does  not  hold  for 
low  density  bottoms. 

We  found  that  structural  variations  in  the  bottom  make  a  larger  contri¬ 
bution  to  intensity  variance  than  topographical  variations,  when  the  aspect 
ratio  is  small,  less  than  about  five.  For  aspect  ratios  more  than  ten, 
topographical  roughness  has  the  greater  effect.  When  comparing  effects  of  the 
different  constituents  of  bottom  topography,  it  was  found  that  curvature 
causes  significantly  greater  intensity  fluctuations  than  either  depth  or  slope 
deviations.  This  is  a  consequence  of  the  tendency  of  curvature  to  focus  and 
defocus  ray  bundles.  Lastly,  small-scale  roughness  is  superimposed  upon  the 
large-scale  facets,  and  we  examined  the  effect  on  intensity  induced  by  this 
feature.  It  is  found  for  our  model  that  only  for  aspect  ratios  less  than  about 
five  does  the  added  roughness  significantly  contribute  to  intensity  fluctua¬ 
tions.  The  small-scale  roughness  tends  to  increase  variance  of  intnesity. 
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"Actual"  and  "Ideal"  rays  for  large-scale  interface  variations. 

Moment  ratio  M^  versus  mean  bottom  density  p2  for  four  L'  combinations: 
R1  =  20,  or  =  a, =  o  ,  =  1,  c  =  1523  m/s,  p,  =  1.025  g/cm3.  Surfaced 
source  and  receiver. 

Moment  ratio  M  versus  mean  bottom  density  p  for  three  values  of 
K  ^ 

aspect  ratio  R' :  =  L'^,  =  L^,  =  0.  Other  parameters  as  in  Fig.  2. 


Moment 


ratio  Mj^  versus  aspect  ratio  R'  for  three  combinations  of  random 


variables:  p2  =  1.8  g/cm3.  Other  parameters  as  in  Fig.  3. 

Same  as  Fig.  4,  except  p2  =  1.6  g/cm.3 

Moment  ratio  M  versus  the  standard  deviations  of  bottom  structure 
K 

o£  and  curvature  o^, :  R'  =  10,  p2  =  1.8  g/cm3.  Other  parameters  as 
in  Fig.  3‘ 

Moment  ratio  Mp  versus  mean  bottom  density  p2  and  aspect  ratio  R' . 
Other  parameters  as  in  Fig.  3. 

Moment  ratio  MR  versus  aspect  ratio  R'  for  three  values  of 
(/:  p2  =  1.8  g/cm3.  Other  parameters  as  in  Fig.  3. 
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